Reconstruction of solar irradiance using the 
Group sunspot number 



o 
o 

(N 



L. Balmaceda, N. A. Krivova and S. K. Solanki 

Max- Planck- Insitut fur Sonnensystemforschung, 
Max- Planck- Str. 2, 37191 Katlenburg-Lindau, Germany 
e-mail: balmaceda@mps.mpg.de 



> 

o 
o 

9 L,: 
6 

CO 



Abstract 

We present a reconstruction of total solar irradiance since 1610 to the present based 
on variations of the surface distribution of the solar magnetic field. The latter is cal- 
culated from the historical record of the Group sunspot number using a simple but 
consistent physical model. Our model successfully reproduces three independent 
data sets: total solar irradiance measurements available since 1978, total photo- 
spheric magnetic flux from 1974 and the open magnetic flux since 1868 (as empiri- 
cally reconstructed from the geomagnetic aa-index) . The model predicts an increase 
in the total solar irradiance since the Maunder Minimum of about 1.3 Wm~ 2 . 
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1 Introduction 



Total solar irradiance changes by about 0.1% between solar activity maximum 
and minimum. Accur ate measurements of this quantity are only available since 
1978 (jFrohlichl . 120061 ) and do not provide information on longer-term secular 
trends. In order to reliably evaluate the Sun's role in recent global climate 
change, longer time series are, however, needed. They can only be assessed with 
the help of suitable models. The most successful models are those attributing 
irradiance variations on tim e-scales longer t han a day to the evolution of th e 
Sun's surface magnetic field (jSolanki et all l2005t Krivova and Solankil . 1200a ) . 
Such models explain more than 90% of all observed ch a nges i n the total and 



spectral irradiance at these tim e scales (IKrivova et all 120031 ; IWenzler et al. 



2004 . l200a . 120061 ; iKrivova et all 120061 ). The continuously evolving distribu- 
tion of the solar magnetic field on the surface is described in these models by 
recourse to magnetograms, which are only available since 1974. For a longer 
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term reconstruction, another proxy of solar magnetic activity has to be em- 
ployed. The available historical proxies of solar activity, such as the Group and 
Zurich sunspot numbers, sunspot, facular or Ca II plage areas mainly describe 
the evolution of the larger magnetic features, such as sunspots or faculae, but 
do not provide any direct information about the weaker features. Therefore, 
whereas the reconstruction of the cyclic component of the irradiance variation 
is typically not a proble m, evaluation of the amoun t of the secular change is 
not straightforward (see ISolanki and Krivoval . 12004 . and references therein). 



Solanki et aD tod l2002h proposed a simple physical mechanism which can 



lead to such a secular trend in the magnetic flux based on the overlap of consec- 
utive activity cycles. Here, we use their model to reconstru ct the magnetic flux 
of th e Sun back to 1610 from the Group sunspot number (IHoyt and Schattenl . 



19981 ) which is employed to reconstruct total solar irradiance for the same 
period. 



2 Approach 



2.1 Photospheric Magnetic Flux 



The basic assumption of our model is that the irradiance variations are caused 
entirely by t he evolution o f the m agnetic features on the solar surface. As in 
the model of Solanki et al.l ( 2002 ). magnetic features on the Sun's surface are 
divided into active regions (AR) and ephemeral regions (ER). The flux emer- 
gence rate in AR, <p act , can be estimated from the Group sunspot number 
since it serves as a good proxy for the fresh flux threading the solar surface. 
The time evolution of the flux emerging in ephemeral regions, <f) e ph, is more 
uncertain, however. Observations suggest that the emergence rate in ER is 
related to that of AR, but the exact shape of this relationship is not yet well 
established. They show that ER associated with the new cycle start emerging 
at the solar surface before the corresponding AR cycle b egins and while mag - 
netic features from the previous cycle are still appearing flHarvevl . [T992l . ll993h . 
Thus, the ER cycle length is extended with respect to that of AR. Therefore 
we prescribe a sine approximation for the shape of the ER cycle, with its length 
being somewhat stretched in time with respect to that of the corresponding 
AR cycle and with its amplitude b eing proportional to the amplitude of the 



AR cycle (see iKrivova et al.l . 120071 . for details). 



Both active and ephemeral regions contribute to the open flux 



'open ) i 



which 



is dragged by the coronal gas and reaches far into the heliosphere. Since this 
flux is mainly unipolar, it decays slowly and lives much longer on the solar 
surface (up to several years) than the AR and ER flux. 
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The extended length of the ER cycles and the long lifetime of the open flux 
lead to an overlap between consecutive cycles such that some background 
magnetic flux is present on the solar surface even at activity minima. The 
amount of this flux changes with time due to variations in the length and 
the amplitude of the magnetic activity cycle. This mechanism thus provides a 
physical explanation for a secular change in the total photospheric magnetic 
flux. 



2.2 Variations of the total solar irradiance 



Following iKrivova et al.l (120031 ); IWenzler et al.l (120051 ). the solar photosphere 
is divided into 5 atmospheric components: the quiet Sun, sunspot umbrae 
and penumbrae, faculae and the network, denoted with subindices q, u, p, f 
and n, respectively. The model consists of two main ingredients: one which is 
temporally invariant and another that introduces a variation with time. 

The time-independent brightness of each component F qiUtP j tn (X), with A being 
the wavelength, is calculated usi ng the ATLAS9 cod e of Kurucz from plane- 



parallel model atmospheres (see lUnruh et al.l . Il999l . for a description of the 
models). Faculae and the network are described by the same model atmo- 
sphere. All fluxes obtained in this way depend only on the wavelength. 

On the other hand, variability in time is due to the changing surface dis- 
tribution of the magnetic components. To describe this, we need to deter- 
mine which part of the solar surface is covered by each component at a given 
time, i.e. the corresponding filling factors, a. In case of sunspots they are di- 
rectly extracted f rom t he sunspot area time series available since 1874 (see 
Balmaceda et al.l . 120051 ) . Before that time, sunspot areas are extrapolated by 
comparing them with the sunspot number. To estimate the filling factors for 
umbrae and penumbrae separatel y, we use the umbral t o penumbral area ratio 
a u /(a„ + a p ) = 0.2 as found by IWenzler et al.l (120061 ). The filling factors of 
the other components are obtained from the reconstructed solar magnetic flux. 
The magnetic flux in faculae, <ftf, can be obtained from (f)t 



-'act 



where 



■'act 



is the flux in AR and 



represent the flux in sunspot umbrae 



and penumbrae, respectively. Finally, the evolution of the network magnetic 
flux, <f> n , which is responsible for the secular change, is given by the sum of 



the flux from ER, (j) ep h, and the open flux, 



■'open- yn 



^eph 



'open- 



The final step is the conversion of the magnetic flux in faculae and the net- 
work into the corres ponding filling facto rs. F or this, we foll o w the conversion 
scheme described by iFligge et al.l (120001 ) and IKrivova et al.l (120031 ). The facu- 
lar and network filling factors increase linearly from at = to 1 at 4> sa t- 
For magnetic flux larger than 4> sa t, the filling factor remains unity. Follow- 
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ing iKrivova et aD (l2003h : IWenzler et aD fl2004j . 120051 . l2006t ) we use the value 



4>satj — 300 G for faculae while for the network we employ (j) sa t,n — 500 G . 
This somewhat enhanced saturation level for the network takes into account 
the fact that a significant amount of weak magnetic flux in the network is lost 
due to the insufficient spatial resolution and relatively high noise level of the 
magnetograms employed. 

Finally, the solar radiative flux at a given wavelength, A, can be obtained by 
combining the fluxes from the 5 components: 



F(X,t) = a q (t)F q (X) + 
+ a u {t)F u {\) + 
+ a p (t)F p (X) + 
+ (a f (t) + a n (t))-F f (X). 



Here, a q (t) — 1 — a u (t) — a p (t) — a n (t) — ct/(t). The total solar irradiance is 
then obtained by integrating F(X,t) over all wavelengths. 



3 Results 



The reconstructed total magnetic flux for individual Carrington rotations is 
compared to data obtained at different observatories, WSO, NSO KP and 
MWO, in Fig. la. The model reproduces both the amplitude and the length 
of the solar cycle in the observed magnetic flux. The evolution of the magnetic 
flux in active and ephemeral regions, as well as of the open and total magnetic 
flux is shown in Fig. l b. The effect of the spatial r esolution on the detection 
of small ER noted by IKrivova and Solankil (120041 ) is taken into account by 
considering the quantity 4> to t = (pact + 0.4 ■(fi ep h + 4 ) open- Since ER cycles overlap, 
the flux in these regions varies only by a factor of 2 over a cycle, which is much 
smaller than the variation in the AR cycles. The ER flux is comparable to that 
of AR during a ctivity maxim a whi le it clearly dominat e s duri ng; minima, in 
agreement with iHarveyl (119931 ) and IKrivova and Solankil (120041 ) . 



Our model also repro duces the secular increa se in the open flux during the 
last century found by lLockwood et al.l (119991 ). This is illustrated in Fig. 2, 
where the modelled open flux is compared to the reconstruction based on the 
aa-index by these authors. 



Figure 3a shows the comparison between the PMOD composite of TSI (jFrohlichl . 



2006 



grey solid line) and the solar irradiance reconstructed from the Group 
sunspot number (dotted line). In order to facilitate comparison, we plot 3- 
month running means. Note, however, that individual dips due to sunspots 
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are also reproduced if daily time series is considered. Both the amplitude and 
the phase of the variations are reasonably reconstructed. Of course, the coarse 
magnetic flux model we use is rather simple to reproduce all details of the 
irradiance variations. For example, the sine approximation of the cycle shape 
does not produce the double-peak structure of cycle 23, as was observed. 

The daily values of the reconstructed solar irradiance since 1610 are shown 
in Fig. 3b. The 11-yr running mean is indicated by the grey solid line. As 
mentioned before, for the period prior to 1874 when no sunspot area mea- 
surements are available, they are obtained by extrapolating sunspot number 
back to 1611. For the period prior to 1753 only monthly values of the Group 
sunspot number can be used. After 1753 daily values are available although 
prior to approximately 1850 there are many gaps in the data. Sampling be- 
comes more regular after 1818. Our model predicts an increase in the solar 
irradiance since the end of the Maunder Minimum (i.e., 1700) till the present 
(average over about 30 years) of about 0.095% or 1.3 Wm~ 2 . 



4 Conclusion 



We have reconstructed total solar irradiance back to 1610. The cyclic varia- 
tion of ER was assumed to be related to the properties of the corresponding 
AR cycle, whose var iation can be estimated from the Group sunspot number 



(jSolanki et al.l . 120021 ). The secular change in the total magnetic flux of the Sun 



and, therefore, in the irradiance is caused by the overlap of the consecutive 
ER cycles. The predicted secular change since 1700 is about 1.3 Win -2 . This 
value lies within the range suggested by oth er recent reconstructions of solar 



irradiance (jFosterl . 120041 ; I Wang et al.l . 120051 ). but is significantly lower than 



the ones obtained in earlier inv e stigat i ons based o n stellar data ranging from 



2 to 16 Wm 2 (e.g., iMendozal . 119971 ; iLeanl . 120001 ) . However, the stellar evi- 



dence for such a change has been rece ntly critized (IHall and Lockwoodl . 120041 ; 



Wright et all 12004 ; iGiampapal . 120051 ) and the magnitude of the increase in 



TSI obtained using these results might have been overestimated. 
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Fig. 1. (a): Total magnetic flux for individual Carrington rotations between 1974 and 
2002. Differen t symbols represe nt data from synoptic charts obtained by different 



observatories (lArge et all |2002j) . The solid line shows the magnetic flux calculated 



from the Group sunspot number, (b): Reconstructed magnetic flux in active regions 
(grey line), ephemeral regions (dotted), as well as the open (dot-dashed) and total 
flux (black line). 
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Fig. 2. The calculated o pen flux (solid line) an d the reconstruction based on the 
geomagnetic aa-index by Lockwood et al. ( 1999I ). 
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Fig. 3. (a): 3-month running mean of the daily sampled TSI for recent cycles: mod- 
elled (dots) and measured (solid line; PMOD composite), (b): Variation of the total 
solar irradiance since 1610. The grey solid line represents the 11-yr running mean 
of this variation. 
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